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ABSTRACT 

High-resolution spectroscopic searches for the starhght reflected from close-in extra- 
solar giant planets have the capability of determining the optical albedo spectra and 
scattering properties of these objects. When combined with radial velocity measure- 
ments they also yield the true mass of the planet. To date, only two such planets have 
been targeted for reflected-light signals, yielding upper limits on the planets' optical 
albedos. Here we examine the prospects for future searches of this kind. We present 
Monte Carlo estimates of prior probability distributions for the orbital velocity ampli- 
tudes and planet / star flux ratios of 6 bright stars known to harbour giant planets in 
orbits with periods of less than 5 days. Using these estimates, we assess the viability 
of these targets for future reflected-light searches using 4m and 8m-class telescopes. 
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1 INTRODUCTION 

Spectral models of close-orbiting giant exoplanets (Seager 
& Sasselov 1998; Marley, Gelino & Stephens 1999; Gouken- 
leuque, Bezard & Lellouch 1999; Sudarsky, Burrows & Pinto 
2000; Barman, Hauschildt & AUard 2001; Sudarsky, Bur- 
rows & Hubeny 2003; Baraffe et al. 2003) show that scat- 
tered starlight dominates over thermal emission at optical 
wavelengths, but the albedo is very sensitive to the depth 
of cloud formation. Sudarsky, Burrows & Pinto (2000) find 
that the relatively low surface gravities of planets such as 
V Andromeda b, HD 75289 b and HD 209458 b may favour 
the formation of relatively bright, high-altitude silicate cloud 
decks. On higher-gravity objects such as r Bootis b, how- 
ever, the same models predict a much deeper cloud deck. 
In this case, much of the optical spectrum is absorbed by 
the resonance lines of the alkali metals. The Na I D lines 
in particular are strongly broadened by collisions with H2 
molecules. 

In reality, however, the atmospheric structure of close- 
orbiting giant planets (hereafter "Pegasi Planets") is ex- 
pected to be more complex, with strong winds and hori- 
zontal temperature variations (Showman & GuiUot 2002), 
and possible temporal variations (Cho et al. 2003). The 
combined atmospheric circulation and temperature varia- 
tions are likely to yield a disequilibrium chemical compo- 
sition. Furthermore, the location and characteristics of any 
cloud decks cannot be predicted to any great precision with 
present theoretical knowledge. The next logical step towards 
understanding the properties and evolution of these objects 
will be to test these emerging models by measuring the albe- 
dos of these Pegasi planets directly. 



A planet of radius Rp orbiting a distance a from a star 
intercepts a fraction {Rp/2af' of the stellar luminosity. The 
proximity of Pegasi planets to their parent stars make them 
excellent prospects for reflected-light searches. At optical 
wavelengths, starlight scattered off the planet's atmosphere 
is expected to dominate over thermal emission. Charbon- 
neau et al. (1999); Collier Cameron et al. (1999) showed 
that the planet/star flux ratio 
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is expected to vary with the star-planet-observer illumina- 
tion angle a, giving a small periodic modulation to the sys- 
tem brightness in the form of the phase function g{a,X). 
Here p(A) represents the geometric albedo of the planet's 
atmosphere that, like the phase function, is wavelength (A) 
dependent. Assuming a grey albedo model, i.e. no wave- 
length dependence, the amplitude of optical flux variability 
for a typical Pegasi planet system is expected to be 



A/ 



< 8.3 X 10 



a/0.045AI7 



(2) 



Future space-based photometry missions are expected to be 
able to detect and measure this modulation with ease. How- 
ever, the orbital inclination (and hence planet's mass) can- 
not be determined directly form the light-curve, since the 
form of the phase function g{a,X) is unknown. 
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Star 


Spectral 
Type 


mv 
(Mags) 


(K) 




Prot 

(Days) 


i; sini 
{km.s~^) 


[Fc/H] 




Age 
(Gyrs) 


M./Mq 




R^/Rq 




V And 


F8V 1 


4.09 


6107±80 


4 


12.2±2 


9.3±0.4 


0.09±0.06 


4 


3.8±1.0 


1.27±0.06 


14 


1.67±0.06 


14 


r Boo 


F7V 1 


4.50 


6360±80 


4 


3.3±0.1 ^ 


14.8±0.3 


0.27±0.08 


4 


l.OitO.6 


1.37±0.05 


14 


1.46±0.06 


14 


51 Peg 


G2V 1 


5.46 


5793±70 


4 


29±7 "5 


2.1±0.7 


0.20±0.07 


4 


4.0±2.5 ■* 


1.10±0.06 


14 


1.20±0.07 


14 


HD 179949 


F8V 2 


6.25 


6115±50 


5 


9±2 8 


6.3±0.9 9 


0.22±0.07 


11 


3.5±2.5 13 


1.24±0.10 


14 


1.24±0.10 


14 


HD 75289 


GOV ^ 


6.35 


6000±50 


3 


16±3 ^ 


4.4itl.O 1" 


0.29±0.08 


12 


5.6±1.0 12 


1.22±0.05 


14 


1.25±0.05 


14 


HD 209458 


F9V 1 


7.65 


6000±50 


8 


16.5±3 1* 


3.7±1.3 


0.00±0.07 


15 


5.7±1.0 16 


l.lOitO.lO 


15 


1.20±0.10 


15 



Table 1. Stellar parameters for bright stars known to harbour Pcgasi planets. ^ Gray, Napier & Winkler (2001), ^ Houk & Smith- 
Moore (1988), 3 Gratton, Focardi & Bandiera (1989), -^Fuhrmann, Pfeiffer & Bernkopf (1998), ^ Gray (1992), Henry et al. (2000a), 
Synchronised rotation assumed from Baliunas et al. (1997); Henry et al. (2000a), ® Mazeh et al. (2000), ^ Groot, Piters & van Paradijs 
(1996), 1° Benz & Mayor (1984), Tinney et al. (2001), Udry et al. (2000), Estimate based on spectral type Average of 
Spectroscopic studies by Fuhrmann, Pfeiffer & Bernkopf (1997); Fuhrmann, Pfeiffer & Bernkopf (1998); Gonzalez (1998); Gonzalez & 
Laws (2000); Gonzalez et al. (2001) Estimated from Tinney et al. (2001); Udry et al. (2000), Cody & Sasselov (2002) i'' Queloz 
et al. (2000), i* Average from Barnes (2001); Queloz et al. (2000) 



Planet 




^orb 


Orbital 


Mp/Mj^p sini 




Rp 1 R Jupiter 


Rp / R Jupiter 




(m.s-i) 


(Days) 


Radius (au) 




(K) 


Hot model 


Cold Model 


V And b 


74.5±2.3 1 


4.6171±0.0001 1 


0.0588±0.0020 


0.716±0.053 


1570 


1.49 


1.32 


T Boo b 


469±5 2 


3.3125±0.0002 i 


0.0483±0.0018 


4.242±0.224 


1680 


1.37 


1.16 


51 Peg b 


56±1 3 


4.2306±0.0005 


0.0528±0.0029 


0.475±0.038 


1330 


1.53 


1.48 


HD 179949 b 


102±3 


3.0930±0.0001 * 


0.0446±0.0036 


0.838±0.109 


1550 


1.40 


1.31 


HD 75289 b 


54=tl 5 


3.5098±0.0007 ^ 


0.0483±0.0020 


0.461±0.028 


1470 


1.58 


1.70 


HD 209458 b 


81±6 6 


3.5239±0.0001 ^ 


0.0470±0.0018 


0.620±0.050 


1460 


1.38 


1.36 



Table 2. System parameters for known Pegasi planets around bright stars. ^ Private Communication - Geoff Marcy, ^ Butler et al. 
(1997), 3 Mayor & Queloz (1995); Marcy et al. (1997), ^ Tinney et al. (2001), ^ Udry et al. (2000), ^ Henry et al. (2000b). The results 
in columns 4 and 5 are calculated from stellar and planetary parameters using the equations detailed in Section 3. Effective planetary 
temperature and radius estimates are given for the estimated age of the system, assuming an edge-on inclination i = 90 (i.e. upper limit). 



2 SPECTROSCOPIC REFLECTED-LIGHT 
SEARCHES 

Starlight reflected from a planet's atmosphere contains 

copies of all the thousands of photospheric stellar absorp- 
tion lines, Doppler-shiftcd by the planet's orbital motion 
and modulated in strength by the phase function g{a,\). 
By detecting and characterising the planetary refiected-light 
signature we observe the planet/star flux ratio as a function 
of orbital phase and wavelength. The information we aim to 
obtain (in order of increasing difficulty) comprises: 

• Kp, the planet's projected orbital velocity. From this we 
learn the orbital inclination i and hence the planet mass Mp, 
since Mp sin i is known form the star's Doppler wobble. 

• eo, the maximum strength of the reflected starlight that 
would be seen if the planet were viewed fully illuminated. 

• p(A), the albedo spectrum, which depends on the compo- 
sition and structure of the planet's atmosphere. 

• p(a,A), the phase function describing the dependence of 
the amount of light reflected toward the observer by the 
star-planet-observer angle a. 



2.1 Current Status 

Several attempts have been made to detect this "faint" 
echo of the starlight, using high- precision spectral subtrac- 
tion and pattern-matching techniques. Charbonneau et al. 
(1999); Collier Cameron et al. (1999) and Leigh et al. (2003) 
searched independently for starlight reflected from r Bootis 
6, establishing an upper limit on the planet/star flux ratio 
that would be seen if the planet were viewed fully illumi- 
nated, with 



eo 



< 1Q~ 



(3) 



for the observed wavelength range of 407 - 649 nm. For an 
assumed planetary radius Rp = 1.2 Rjup, the Leigh et al. 
(2003) 99.9% confidence levels (Fig. 3a) suggest a geometric 
albedo p < 0.39 at the most probable orbital inclination 
i ~ 40° , assuming a grey albedo model. 

Collier Cameron et al. (2002) used the same technique 
on V Andromeda h, establishing an upper radius limit Rp < 
1.22Rjup, assuming a grey geometric albedo of p = 0.5. 
(Fig. 3b). Both Collier Cameron et al. (2002) and Leigh 
et al. (2003) find possible detections above the 2a signifi- 
cance level, which if confirmed would yield masses of 0.74 
and 7.28 Mjup respectively for v And b and r Bootis 6. How- 
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ever, there is a ~ 4 to 8% probability that the detected fea- 
tures, shown at Figure 3, are a consequence of non-gaussian 
noise, and as a result neither are confident enough to claim 
a genuine detection. 



3 SYSTEM PARAMETERS 

The known stellax parameters for 6 bright stars known (from 
Doppler wobble studies) to harbour Pegasi planets are set 
out in Table 1. Fits to the radial velocity data have pro- 
vided us with orbital periods for the planets, together with 
the stellar reflex velocity of each star about the com- 
mon ccntre-of-mass of the system. Both parameters are set 
out in Table 2 alongside calculated values for the orbital 
radius of the system (from Kepler's Laws) and the mini- 
mum planet mass. In addition, we provide theoretical es- 
timates for the effective temperature and upper limit on 
the radius of the planet (for edge-on orbits), assuming both 
"Ifot" and "Cold/Dissipative" atmospheric structure mod- 
els, as described in Section 3.2. 



3.1 Planetary Mass Estimates 

Among the known Pegasi planets there are several promising 
candidates for future spectroscopic reflected-light searches. 
The main selection criteria are that the host star must be 
bright, and that the planet/star flux ratio near superior con- 
junction must be high. Given spectroscopic and/or photo- 
metric estimates of the stellar radius 7?,, the projected ro- 
tation speed vsini and rotation period Prot, we estimate 



27ri?* 



(4) 



This makes the assumption that the stellar rotation axis and 
the orbital plane are close to perpendicular, as was deter- 
mined for HD 209458 by Queloz et al. (2000). The mass ratio 
q — Mp/M-M follows from the observed orbital period Porb 
and stellar reflex velocity if* : 



K.P„, 



1 + q 2'Ka sin i 



(5) 



This then provides us with the planet's projected orbital 
velocity amplitude {Kp = Kt/q) in addition to the planet 
mass {Mp = qMt). 



3.2 Planetary Radius Estimates 

Assuming a known mass, the radius of a gaseous planet is 
governed by its Kelvin-Helmholtz contraction, and the fol- 
lowing factors, in order of significance: (1) its composition; 
(2) its atmospheric temperature; (3) its age; and (4) possible 
additional energy sources. Uncertainties also stem from our 
limited understanding of the underlying physics: equations 
of state, opacities - see Guillot (1999). 

The measurement of the radius of HD 209458 b (Char- 
bonncau ct al. 2000; Henry et al. 2000b; Brown ct al. 2001) 
showed that the planet is essentially made of hydrogen and 
hehum (Guillot et al. 1996; Burrows et al. 2000), thus ad- 
dressing the first source of uncertainty concerning the radii 



of Pegasi planets. However, it was later shown that the ra- 
dius of HD 209458 b is too large to be explained by stan- 
dard evolution models using realistic atmospheric bound- 
ary conditions (Bodenheimer, Lin & Mardling 2001; Guillot 
& Showman 2002; Baraffc ct al. 2003). Aflhough Burrows, 
Sudarsky & Hubbard (2003) question the need to invoke 
additional atmospheric energy sources, Bodenheimer, Lin & 
Mardhng (2001) proposed that HD 209458 b may be kept in- 
flated by the dissipation of tidal energy through the circular- 
izatiori of its orbit. Due to the very small eccentricity of the 
orbit (consistent with 0.0), Showman & Guillot (2002) pro- 
posed two alternative explanations: either the atmosphere 
is relatively cold, as predicted by radiative transfer calcu- 
lations, and kinetic energy generated in the atmosphere is 
dissipated relatively deep in the interior (e.g. by tides due 
to a slightly asynchronous rotation of the interior), or the 
atmosphere is hotter than calculated, for example because 
shear instabilities force it into a quasi-adiabatic state. 

We therefore calculate the radii of other Pegasi planets on 
the following basis: 

• All Pegasi planets have approximately the same com- 
position (they are mostly made of hydrogen and helium) . 

• Their atmosphere is either "hot" or "cold" , the temper- 
ature at a given level (10 bars) being calculated as a func- 
tion of the effective temperature and gravity, as described in 
Guillot & Showman (2002) (see also Burrows et al. (1997)), 
assuming a bond albedo equal to 0.4. 

• In the "cold" case, an additional energy flux is deposited 
at the centre of the planet. This flux is set equal to 8 x 
10~^ times the absorbed stellar luminosity (i.e. between 1.2 
and 3.1 X lO^^ergs"^). This corresponds to the energy flux 
necessary to reproduce the radius of HD 209458 b with the 
"cold" boundary conditions (Guillot & Showman 2002). 

Of course, a number of unknowns remain. Most impor- 
tantly, both the composition, bond albedo and the energy 
flux dissipated in the planet are likely to be complex func- 
tions of the planetary mass, orbital distance and even stellar 
type. However, this is likely to be a second order effect be- 
cause our radii calculations take advantage of the known 
radius of HD 209458 b. It is interesting to see that even 
though the parameters of the models have been set by the 
constraint on HD 209458 b, the "hot" and "cold" scenarios 
shown in Figure 1 yield signiflcantly different radii for the 
other planets. This is due to the fact that energy dissipa- 
tion (in the "cold" case) has a greater impact on planets of 
smaller mass, and that the "hot" boundary condition yields 
a planetary radius that is more sensitive to the orbital dis- 
tance. 



3.3 Flux Ratio Estimates 

For a Lambert-sphere phase function g{a) we can compute 
the maximum expected planet /star flux ratio near to supe- 
rior conjunction, 

' Rp ^ ^ 



~a) 



(6) 



We adopt the Lambert-sphere phase function because it pro- 
vides a good approximation to the phase functions of gas gi- 
ants within our own solar system (Pollack et al. 1986; Char- 
bonneau et al. 1999). Although we expect the albedo spectra 
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Figure 1. Evolutionary radius models for the planets surrounding v Andromeda, t Bootis, 51 Pegasi, HD 179949, HD 75289 and HD 
209458. The evolutionary tracks show both the "Hot" and "Cold/Dissipativc" models for orbital inclinations of sin? = 0.5 and sini = 1.0 
(Edge-on). Being a known transiting planet, HD 209458 b is shown solely for sini = 1.0. 



of Pegasi planets to exhibit a strong wavelength dependence, 

here we adopt a simple grey geometric albedo p = 0.4 as a 
plausible comparison to the case of a high, reflective sil- 
icate cloud deck with little overlying absorption, as in the 
Class V models of Sudarsky, Burrows & Pinto (2000). Other 
models in which cloud layers are deeper, or even absent, 
produce lower optical albedos; thus the flux ratios we com- 
pute here should be treated as upper limits on the plausible 
planet/star flux ratio. 



4 METHOD LIMITATIONS 

Spectroscopic reflected-Ught searches should provide us with 

a unique opportunity to achieve the direct detection of a 
Pegasi planet, however there exist a rmmber of potential 
difficulties that require a cautious approach to any analysis. 

(i) At the flux ratios we are working {fp/f* < 10"*), 
we inevitably see the appearence of systematic non-gaussian 
noise within data sets, noise that requires to be carefully 
identified and corrected for in the subsequent analysis. 

(ii) Any predictions from our searches require assump- 
tions to be made about the nature of the planetary system, 
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V Andromedo b t Bootis b 51 PegasI b 




50 100 150 200 50 100 150 200 50 100 150 



Kp km.s ' Kp km.s ^ Kp km.s 



HD 179949 b HD 75289 b HD 209458 b 




50 100 150 200 50 100 150 200 50 100 150 



K km.s ' K km.s ' K km.s ' 

p p p 

Figure 2. Prior probability density maps for tlic plaiict/star flux ratio near superior eonjunctioii, as a function of the planet's projected 
orbital velocity Kp, for an assumed geometric albedo p = 0.4. The maps are calculated for an evolutionary radius model that exhibits a 
"Hot" atmospheric structure, although maximum likelihood values for both the "Hot" and "Cold" scenarios are shown at Table 3. Trials 
in which the orbital inclination was high enough to cross the observer's line of sight to the stellar disc were rejected for all cases bar HD 
209458 b, since none of these other systems has been seen to exhibit transits (Henry et al. 2000a; Tinney et al. 2001). For the transiting 
planet, its known inclination and mass act to constrain the projected radial velocity amplitude. 

In the case of rBootis, an additional constraint was applied. Since the rotation of the primary appears to have become synchronised with 
the planet's orbit (Baliunas et al. 1997; Collier Cameron et al. 1999; Henry et al. 2000a), the planet must have been massive enough to 
have synchronised the star's spin within its own main sequence lifetime. For each trial we computed both the main sequence lifetime and 
the synchronisation timescale. Trials in which the latter exceeded the former were rejected. This, in conjunction with the well determined 
vsini value, leads to a fairly tight constraint on the system inclination and hence Kp. 



Planet 


Kp 
{km.s~^) 


Inclination 
(Degrees) 


Mp/Mj^p 


ifp/ }*)max 

Hot model 


(.fp/ f*)max 

Cold Model 


V And b 


133 


78 


0.760 


4.88 X 10"^ 


4.53 X 10-5 


r Boo b 


95 


37 


6.813 


4.18 X 10-5 


3.22 X 10-5 


51 Peg b 


117 


63 


0.523 


4.81 X 10-5 


4.46 X 10-5 


HD 179949 b 


135 


60 


0.945 


6.44 X 10-5 


5.34 X 10-5 


HD 75289 b 


141 


75 


0.498 


7.38 X 10-5 


9.17 X 10-5 


HD 209458 b 


144 


87 


0.620 


7.18 X 10-5 


7.17 X 10-5 



Table 3. Maximum likelihood statistics for the parameter distributions returned by the Monte Carlo trials, given what we already know 
about the system from observation. In adopting a grey albedo of p = 0.4 for our trials, similar to that of the Class V models of Sudarsky, 
Burrows & Pinto (2000). In reality, we would expect the albedo spectra of Pegasi planets to exhibit a strong wavelength dependence. 
However, we believe the use of a grey albedo spectra provides us with a realistic upper limit to the expected planet/star flux ratio. 
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Upper limits Upper limits 




Kp km/sec Kp km/sec 

Figure 3. Results from rcflcctcd-light studies of (a) r Bootis (Leigh ct al. 2003) and (b) v Andromeda (Collier Cameron ct al. 2002), 
shown as relative probability maps of model parameters Kp and log(e()) = ^ogp(Rp/ay^. The results are derived from WHT/UES 
observations of r Bootis and v Andromeda, assuming a grey albedo speetrum. The greyseale denotes the probability relative to the 
best-fit model, increasing from for white to 1 for blaek. The eontours show the confidenee levels at which candidate detections can 
be ruled out as being caused by spurious alignments of non-Gaussian noise features. Prom top to bottom, they show the 99.9%, 99.0%, 
95.4% and 68.3% confidence limits. The dark features in each probability map represent candidate planetary detections from the data 
but occur at too low a confidence level to guarantee a genuine detection. We note, however, that in both cases the main features appear 
close to the Kp values predicted at Table 3 by Monte Carlo simulations and shown at Figure 2. 



Target mv [fp/ f*)max 4m-class 8m-class 

(Mags) Hot model telescopes telescopes 



V And 


4.09 


4.88 


X 


10" 


-5 


20.0 


hrs 


5.0 hrs 


T Boo 


4.50 


4.18 


X 


10" 


-5 


39.8 


hrs 


9.9 hrs 


51 Peg 


5.46 


4.81 


X 


10" 


-5 


72.7 


hrs 


18.2 hrs 


HD 179949 


6.25 


6.44 


X 


10" 


-5 


84.0 


hrs 


21.0 hrs 


HD 75289 


6.35 


7.38 


X 


10- 


-5 


70.1 


hrs 


17.5 hrs 


HD 209458 


7.65 


7.18 


X 


10" 


-5 


245 


hrs 


61.3 hrs 



Table 4. Exposure time predictor to compare benefits of observing different star and telescope class combinations. We take 4m-class 
observations on v Andromeda as our standard, with an exposure time benchmark of 20.0 hrs required to reach the predicted {fp/ f*)m.ax 
signal levels. In deriving these estimates we assume that all other variables are equal, such as overall instrument efficiency, observing 
strategies and weather conditions. 



such as the phase function g{a,\), the planetary radius Rp 
or albedo spectra p{X), and that the orbital and stellar equa- 
torial planes are co-aligned. To some extent observational 
results will therefore act to narrow the range of parameter 
space occupied by these Pegasi planets, although stronger 
detections will allow us to draw more definite conclusions. 

(iii) The "Doppler Wobble" method of detection con- 
firms the existence of a planet but cannot tell us its true 
orbital inclination. A spectroscopic search for the reflected- 
light component of a close orbiting planet has the potential 
to identify the projected orbital velocity Kp and hence the 
orbital tilt. However, with such a scope of possible Kp val- 
ues, it is important to estimate where we might expect to 
detect the planet, given the known parameters of the system. 
Such estimates, as detailed in Section 5, provide us with a 
quantifiable method by which we can assess the plausibility 
of any candidate detection. 



5 PROBABILITY DISTRIBUTIONS 

We derive the prior probability distributions for the observ- 
able quantities Kp and {fp/ f*)max by drawing randomly 

from the uncertainty distributions (assumed gaussian) for 
the observed quantities AT*, M«, Prot, Age and wsini. 
The values of these quantities, for 6 of the brightest stars 
known to harbour Pegasi planets with orbital periods of less 
than 5 days, are given in Tables 1 and 2. In order to pro- 
vide a radius estimate for each planet Rp, we use our sini 
estimate to conduct a log Rp — log Mp linear interpolation 
(for a selected Age) between the sini = 1.0 and sini = 0.5 
evolutionary models. This is done for both the "Hot" and 
"Cold" evolutionary scenarios. In each case we performed 
1,000,000 random Monte Carlo trials to compute Kp and 
^og{fp/ f,)m,ax- The resulting joint probability distributions 
for the planet/star flux ratio and projected orbital velocity 
amplitude are presented in Figure 2. 
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6 DISCUSSION 

The maximum observable planet/star flux ratio for a planet 
depends strongly on the inclination of its orbit to our line 
of sight, and on the distance of the planet from the star. 
The statistical analyses presented here provide estimates of 
the planet/star flux ratio, derived from the best data cur- 
rently available on 6 bright stars harbouring planets with or- 
bital periods of less than 5 days. These estimates are model- 
dependent, in that we have made theoretical assumptions 
about each planet's radius, the co-alignment of each orbital 
plane to the stellar rotation axis and used a Lambert-sphoro 
phase function to describe the atmospheric scattering prop- 
erties. In addition, we have adopted a grey geometric albedo 
(p = 0.4) that corresponds roughly to the reflective, high- 
level silicate cloud decks predicted by the Class V spectral 
models of Sudarsky, Burrows & Pinto (2000) . 

Recent spectroscopic seaxches for relected-light signa- 
tures from r Bootis b (Leigh et al. 2003) and v Andromeda 
b (Collier Cameron et al. 2002), shown at Figure 3, have al- 
ready reached detection limits comparable to the signal lev- 
els predicted here, using echelle spectrographs on 4m-class 
telescopes. We find the total number of photons that must 
be collected scales with the square of the planet /star flux 
ratio, whilst the length of time required to collect this num- 
ber of photons with a given telescope scales with the flux 
received from the star. If we adopt 4m-class observations 
of V Andromeda as our benchmark (20.0 hrs) to reach the 
ifp/ f*)max signal levels predicted, then Table 4 details the 
relative exposure times for each planetary system, given ac- 
cess to both 4m and 8m-class telescope facilities. Given that 
it has been possible to probe to these deep signal levels for 
V And in only a few nights on a 4m-class telescope, the 
remaining 5 targets are clearly viable and compelling tar- 
gets for future studios with existing high-throughput spec- 
trographs on 8m-class telescopes. Future instruments that 
are also likely to be useful in this respect include the high- 
efficiency fibre-fed echelle spectrograph ESPADONS, which 
is to be commissioned on CFHT during 2003. 
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